Surface Quantized Anomalous Hall Current and Magneto-Electric Effect 
in Magnetically Disordered Topological Insulators 
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We study theoretically the role of quenched magnetic disorder at the surface of a topological 
insulator by numerical simulation and scaling analysis. It is found that all the surface states are 
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localized while the transverse conductivity is quantized to be as long as the Fermi energy is 
within the bulk gap. This greatly facilitates the realization of the topological magnetoelectric effect 
proposed by Qi et al. (Phys. Rev. B78, 195424 (2008)) with the surface magnetization direction 
being controlled by the simultaneous application of magnetic and electric fields . 

PACS numbers: 73.43.-f,75.70.Kw,85.75.-d,85.70.Kh,85.75.-d 



The bulk-surface correspondence has an essential role 
in a large variety of phenomena in condensed mat- 
ter physics, such as ferroelectricity, diamagnetism, the 
Meissner effect, and the quantum Hall effect. The topo- 
logical magnetoelectric (ME) effect is a novel manifes- 
tation of the bulk-surface correspondence in which the 
bulk magnetization is generated by a circulating quan- 
tized Hall current flowing at the surface of topologically 
nontrivial insulators, called topological insulators (TIs)EJ. 

The electromagnetic response of a three-dimensional 
(3D) TI can be depf^ribed by the Lagrangian for the axion 
electrodynamics ,L 
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where E and B are the electromagnetic fields, e and 
fi are the dielectric constant and magnetic permeabil- 
ity, and a — /he is the fine structure constant. The 
first term is the conventional Maxwell term. The second 
term, called the 6'-termEl, characterizes the topological 
nature of three-dimensional insulators; = or tt (mod 
2-k) corresponding to ordinary insulators and topologi- 
cal insulators, respectively. The 9 term is written as the 
total divegence, and hence can be transformed into the 
surface Chern-Simons term in a system terminated by a 
boundary. When a topological insulator is surrounded 
by a ferromagnetic layer, time-reversal (T) symmetry is 
broken, and the value of 9 (more precisely V^) can be 
determined as discussed in Refs. ||,^. A ciwajlating Hall 
current, induced by an applied electric field ip on the sur- 
face, is the source of a bulk magnetizationa. These fea- 
tures arc distinct from the conventional ME cffectEl which 
is a long-term issue in the field of multiferroics. 

Experimentally, however, there are several difficulties 
to realize this topological ME effect, (a) First, it is re- 
quired to get rid of the bulk carriers, (b) Second, one 
needs to attach the insulating ferromagnetic layer with 
the magnetization normal to the surface all pointing out 
or in. (c) Lastly, the Fermi energy must be tuned accu- 
rately within the small gap of the surface Dirac fermion 



opened by the exchange interaction. Otherwise, the de- 
scription by eq.(|l|) is not justified, (a) is progressively 
realizedQ, while (b) and (c) seem still very difficult at the 
moment, even though a recent experimental work shows 
that the gapless surface Dirac states of the pristine topo- 
logical insulator Bi2Se3 become gapped upon introducing 
magnetic impurities (Mn and FCj) into the crystal and also 
Mn dopants lead carrier dopingj. Therefore, the topolog- 
ical ME appears to be not practical c»fjarj-4hough several 
theoretical proposals have been madajollj. 

In this paper, we study the effects of the quenched 
magnetic impurities or disorder on the surface of TL In 
sharp contrast to conventional quantum Hall systems, all 
the surface states are localized while the Hall conductiv- 
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ity is quantized to be as long as the Fermi energy 
is within the bulk gap. This resolves the problem (c). 
Consequently the generated magnetization is robust over 
randomness and universal. Furthermore, it is shown that 
this also resolves (b) with the simultaneous application 
of magnetic and electric fields parallel or antiparallel to 
each other. By this method, the surface magnetization 
can be controlled by the hulk energy, and hence can easily 



Hi„d = (47i/c)yH (b) 




(c) 



FIG. 1: (color online) (a) Illustration of a magnetically doped 
topological insulator with a cylindrical geometry, (b) Surface 
massive Dirac dispersion. The Hall conductivity axy is quan- 
tized when the Fermi level lies in the surface gap. (Jxy deviates 
from e^/2/i when the Fermi level is out of the surface gap. 
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overcome the magnetic anisotropy and Zeeman splitting 
at the surface. 

We study the surface transport of a 3D TI basing 
on the 2D massive Dirac model with various types of 
disorder. The interaction between magnetic dopants 
and surface electrons is described by the exchange 



Hamiltonian 



He 



Ri). The uniform part of z-component of local spins geft- 
erates a mass gap m = JnimpS'z in the surface spectrumnl, 
where cr = ((Tj^, <y2,cr^) are Pauli matrices that act on the 
electron spin degrees of freedom, Si and Ri are the lo- 
cal spin operator and the position operator of ith mag- 



netic dopant, 



iViinp and rii^ 



are the total number and 



the mean sheet density of local spins. The surface |£)irac 
modes can be described by the Dirac Hamiltonianu 



^ofrac = -iflVFZ X CT • V 



(2) 



On the other hand, inhomogeneous part of the local 
magnetization gives disorder term: V — X]|i=o ''^a'^a'(^)' 
where Vi and V2 have a role of random vector potential, 
and V3 has a role of random mass potential. The scalar 
potential Vo(t') also could be introduced by impurities or 
vacancies. Here ctq is a 2 x 2 unit matrix. 

In the clean limit, the Hall conductivity is quantized as 
f^xy = ~sgn(m)|jj when the l^rmiJevel is in the middle 
of the surface gap < |TO|)0iB0. The fact that a^y 
is quantized eveij in the limit of ?n — > is referred to 
parity anomalyu3. When \E\ > |m|, on the other hand, 
axy deviates from the quantized value, Sjod vanishes at 
|i?|/|m| — > 00 as sketched in the Fig.l(c).Ej 

We evaluate the diagonal and the Hall conductivity 
of randmn massive Dirac Hamiltonian with the Kubo 
formulalia 
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c) Uxx l^lh] 



where a, 6 = cc or y, is the current operator, and |n) 
denotes an eigenstatc with its eigenvalue E„ of Hamil- 
tonian Hp^jjf. + V. We work in the momentum space 
by introdu cing , a hard cutoff at a sufficiently large mo- 
mentum A.E3Eil Eigenstates and eigenvalues are obtained 
by numerically diagonalizing the Dirac Hamiltonian with 
disorder terms in the momentum space. Random averag- 
ing is taken over typically 1000-10000 disorder configura- 
tions. We use the Gaussian model for disorder potentials 
which obey {V^{q)VAq')) = S^,g^exp{-q^dy2)S{q+q') 
with ^,1^ = 0,1,2,3. We sejt, go = 5i = ff2 = <?3 as as- 
sumed in realistic situations and the disorder strength 
so that the disorder broadening energy is the order of the 
surface gap m as seen in the density of states shown in 
Fig. 2(a). 

The diagonal Gxx and the Hall conductivity Gxy are 
shown in Fig. 2 for three different system sizes as a func- 
tion of E /m^ E being the Fermi energy. The sign of a^y 
is determined solely by the sign of the mass, no matter 
whether the Fermi level resides in the electron- or hole- 
region. The size dependence indicates that a^x decreases 
while Uxy increases as the size L increases. The linear size 
of the system is characterized by the momentum cutoff 
A or the mass gap m. Three system sizes shown in Fig. 2 
are L ~ 75A^^ = Ahvp/m, 113A~^ = Ghvp/m, and 
176A-1 = 9hvp/m. 

To see the tendency of axy and axx in the L 00 limit, 
we plot (cTxa;, fxy) by changing the system size L in Fig. 
3. Although the system size cannot be widely changed 
because of the computational limitation. Fig. 3 indicates 
that, with increasing system size, {crxxjO'xy) fiows and 
approaches to two fixed points (0, ±1/2) in units of e'^/h. 
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FIG. 2: (color online) (a) Density of states, (b) Hall conduc- 
tivity, and (c) diagonal conductivity of massive Dirac fermions 
are shown as a function of the Fermi energy divided by the 
mass gap E/m. Lines are guides for the eyes. System sizes 
are L = Ahvp /m,6hvF /m, and Qhyp/m 



FIG. 3: (color online) Scaling flow of axx and axy as increas- 
ing system size, indicating {axx, o'xy) — (0, ±1/2), in units of 
e^//i, at L — > 00. Curves are guides for the eyes. Inset: The 
surface Hall current is related to the bulk orbital magnetiza- 
tion induced by an electric field. 
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The scahng law of Uxx and a^y has itufW originally studied 
in the integer quantum Hall effecttJO with correspond- 
ing fixed points (0,ri,) in units of e^//i, n being an inte- 
ger. In the present case with single-flavor Dirac fermions, 
a field theoretical studycJ shows that the values of Uxy 
shifts by 1/2 from the conventional ones, consistent with 
Fig. 3. The 1/2-shift results in a marked point, that is, 
^xy = is unstable. This means that even if the surface 
mass gap is, in principle, infinitely small, \oxy \ increases 
and approaches to jih. This is a generalization of par- 
ity anomaly to the case with disorder. This is in contrast 
to the multi-Dirac-cone 2D lattice systems where mixing 
between different Dpac-cones wipes out the Hall conduc- 
tivity quantization.Es 

When the system size exceeds the localization length, 
the localization effect becomes important and a^y starts 
to flow toward the two fixed points as shown above. The 
maximum system size in above analysis corresponds to 
^ 0.2 [/^m] (with identifying A with the largest wave num- 
ber in the surface spectrum ~ 0.2[A~^]) which is very 
small compared to the realistic sample size and also to 
the coherence length. We expect that in the experimen- 
tal situations the quantization of the Hall conductivity 
should be more prominent. lur-conventional quantum 
Hall systems, recent experimentsc3 explored the temper- 
ature driven fiow diagram of g^x and Oxy in a large- 
temperature range from 4 K down to 40 mK. Similar scal- 
ing properties, but different fixed points, are expected on 
a magnetically doped surface of a topological insulator. 

Wc note that even for states out of the original surface 
gap, which were metallic in the clean limit, {pxx^<^xy") 
scales to the fixed points (0, ±1/2). This means the fact 
that all the surface states are localized, while a trans- 
verse current flows. It has been known, however, that 
in two-dimensional electron gas systems at least one ex- 
tended state below the Fermi level is required to real- 
ize the quantum Hall effeclEj, otherwise Oxy vanishes. 
Since all states of disordered massive Dirac Hamiltonian 
are localized in zero magnetic field, a finite quantized 
Hall transport sounds enigmatic: what carries the elec- 
tric charges? In the following, we address this problem 
basing on the bulk-surface correspondence. 

Since all the surface states are localized by magnetic 
disorder, the description by the 0-term in eq.(|^) is justi- 
fied, which corresponds to the magnetoelectric effect: 



a 



(4) 



These results have been reproduced from the ex- 
plicit WaHTiier state representation of a topological 
insulator .LlrEa Here Q is only well defined as a bulk prop- 
erty modulo 27r. The integer multiple of 27r can be speci- 
fied once we specify a particular way to make the bound- 
ary with broken T-symmetry. More appropriately, the 
current and charge densities are given by 



jf = cV X Me + dPs/dt, p = -V 

consistent with the half-integer qjt|^ 
face states, under the conditio 



Pb, (5) 

Hall sur- 
d9/dz = 



axy{'2TTh/e'^)5{z—R), where z = Ris at the surface shown 
in Fig. 1(a) and the inset in Fig. 3. Since bulk moments 
Eqs.(^ arc protected by the bulk gap (when it remains 
intact) and insensitive to disorder on the surface, Eqs.(|^) 
can survive even though all surface states are localized. 
With above results, the existence of the surface quantum 
Hall states, generalizes the axion electrodynamics effects 
to disordered systems. Precisely, the energy region for 
the topological ME effect, which was originally limited 
in the surface gap, is enlarged to the bulk gap by local- 
ization effect of surface states. This will greatly facilitate 
the experimental realization. 

Even though the 2D surface Dirac Hamiltonian omits 
such the bulk properties by the construction, non- 
vanishing half-quantized Hall conductivity is represented 
as parity anomaly. A similar situation can be seen in the 
context of chiral anomaly. For the one-dimensional edge 
channel of the quantum Hall system, the chiral anomaly 
represents the current fiow from or to 2D bulk region of 
the sample.Ej 

We note that the topological ME effect requires a fi- 
nite surface gap with an entirely definite sign of the mass, 
otherwise there exists a massless one-dimensional channel 
at the,domain boundary separating different signs of the 
mass.Eia In general, however, magnetic moments of impu- 
rities could form a domain structure when temperaiuce 
decreases down below the transition temperaturcLiO. 
Simple situations are sketched in Fig. 4. There are two 
(many in general) domain regions where the surface mag- 
netizations point alternative directions, one S± > (neg- 
ative mass) and the other < (positive mass). 

In the remaining part of the paper, we propose a practi- 
cal way to control the domain structure using the massive 
Dirac surface states. First we assume that a magnetic 
field is applied in +z direction as shown in Fig. 4(a). For 
simplicity, let us start with the pure case without disor- 
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FIG. 4: (color online) Illustration of domain structures and 
electromagnetic control of them. Domain wall is charged in 
the case of (b). 
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der. On the surfaces perpendicular to the magnetic field, 
the Landau level structure is formed. When the Fermi 
energy resides in the mass gap {\E\ < \m\), the electric 
charge density p = ±{e^ /2h)B is generated on the top 
and bottom. The sign of the charge depends on the sign 
of the mass and also the direction of the magnetic field. 
By applying an electric field in z direction, the degener- 
acy between the right and left regions is lifted, because 
of the electric polarization energy. In the positive mass 
region there is an energy gain, while loss in the negative 
region, as estimated as 

~ -10i°(S[V/cm])(B[T])(L[cm])^ [eV], (6) 

where B — H + AnMi„^p includes the contribution from 
the bulk magnetization of impurities Afimp, and L is the 
linear size of the domain region. For this energy gain, 
surface electrons in the right region are transferred from 
top to bottom. Tliis process requires larger Ue than the 
anisotropic energyllj C/aniso/i^ 5 x 10^° [eV/cm^] and 
the Zeeman energy C/zeema„/i' - 109(S[T]) [eV/cm^] 
to fiip the surface magnetization (sign change of the 
mass) in the right region. Typical strengths of electric 
E ^ 10^[V/cm] and magnetic fields B ^ 1 [T] are enough 
to flip and rearrange the surface magnetization. In the 
case of strong disorder so that the Landau level spacing is 
dominated by disorder broadening energy, the extended 
states originally located at the center of each Landau 
level levitate to high (low) energy regime if its energy in 
the clean limit was positive (negative). Ej Eventually, only 
localized states remain in the spectrum where axy scales 
to ±e^/2/i as shown above. In this case, the induced 
charge density on the surface is p ~ cr-xyB independent 



of the Fermi level as long as it is within the bulk gap. 
Consequently, as domain wall moves, the positive mass 
region dominates over the negative mass region, and the 
uniform surface gap with a definite sign is attained, where 
the Hall current can fiow circularly. Once the uniform 
mass rearrangement is formed, the external fields E and 
H can be gradually turned off. When an electric and 
magnetic fields point in x direction (see Fig. 4(b)), a sim- 
ilar effect occurs. The surfaces with generated electric 
charge are different but the direction in which domain 
walls move is same as the case of Fig. 4(a). 

In this work, we studied the effects of quenched mag- 
netic disorder on the surface of a topological insulator, 
basing the surface Dirac model. The scaling analysis 
indicates that all surface states are localized, while the 
Hall conductivity approaches to the quantized value, even 
when the mass gap is smeared by disorder broadening. 
Consequently the plateau width at large system size and 
at low temperature is enhanced from that in the clean 
limit. This helps the experimental realization of the 
topological ME effect in a magnetically doped TI. We 
also proposed the way to control domain structures of 
the surface magnetization by the simultaneous applica- 
tion of electric and magnetic fields. The effective theory 
of TI including the magnetic impurities is an interesting 
and important issue, which is left as a future subject. 
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The massless case has been studied in detail in Ref.^. In 
the massive case, the critical state originally located at 
the center of A'^ = Laudau level levitates to the higher or 
lower energy region, depending on the sign of the mass. 



